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Beschreibung 

Verfahren zura Bearbeiten von Magnetresonanzrohdaten, Magnet- 
resonanztomographiegerat und Magnetresonanzspektroskopiegerat 

Die Erfindung betrifft ein Verfahren zur Nachverarbeitung von 
Magnetresonanzrohdaten. Aufierdem betrifft die Erfindung ein 
Magnetresonanztomographiegerat oder ein Magnetresonanzspekt- 
roskopiegerat, das jeweils an ein solches Verfahren zur Nach- 
verarbeitung von Magnetresonanzrohdaten angepasst ist. 

Die Bearbeitung von Magnetresonanzrohdaten, die sowohl von 
Magnetresonanzspektroskopiegeraten als auch von Magnetreso- 
nanztomographiegeraten stamen konnen, hat zum Ziel, die me- 
dizinisch relevante Information aus den Magnetresonanzrohda- 
ten bestmoglichst zu extrahieren. Im Fall der Magnetresonanz- 
tomographie ist demnach eine Aufnahme mit einer hohen Detail- 
auflosung und niedrigem Rauschen das Ziel, also eine Aufnahme 
mit hoher Kantenscharfe und hohem Signal zu Rauschverhaltnis 
(SNR Signal to Noise Ratio) . Im Fall der Magnetresonanzspekt- 
roskopie soil aus den Rohdaten ein gut aufgelostes Linien- 
spektrum mit hohem SNR gewonnen werden. Die Qualitat in der 
Bearbeitung von Magnetresonanzrohdaten liegt also zum Einen 
in der Kantenscharfe und Bildscharfe, d.h. in scharfen Auf- 
nahmen mit hohem Kontrast, und zum Anderen in der Erzeugung 
von einem hohen SNR, so dass das Wesentliche der Aufnahme 
bzw. des Spektrums hervorgehoben wird. 

Bei den verschiedenen Ansatzen zum Erreichen dieses Ziels 
kann meist eine gleichzeitige Erhohung der Bildscharfe und 
des SNR nicht erreicht werden. Dieser Nachteil wirkt sich di- 
rekt auf die Qualitat der Magnetresonanztomographie und - 
spektroskopie aus. Urn die Detailauf losung zu verbessern, 
iiberlagert man beispielsweise hochf requentes, "blaues" Rau- 
schen einer Aufnahme. Dies erzeugt eine subjektive Erhohung 
der Detailauf losung. Eine ebenfalls subjektive Verbesserung 
der Auflosung bewirkt auf ahnliche Weise die Anwendung von 
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Hochpassf iltern auf die Bilddaten. Allerdings verstarkt auch 
dies das Bildrauschen . 

Um das SNR im Fall von Magnetresonanzauf nahmen zu verbessern, 
werden bekannte Filter (wie beispielsweise im Frequenzraum 
der Harming- , der Fermi- oder der Kosinusf ilter ) auf die 
Magnetresonanzrohdaten angewandt . Die Wirkungsweise solcher 
Filter ist aus der Literatur bekannt, z.B.: F.J. Harris, 
Proc. IEEE, Vol. 66, No. 1 (1978) : „On the Use of Windows 
for Harmonic Analysis with the Discrete Fourier Trans- 
form" . Auch im Ortsraum konnen Nachbearbeitungsf ilter ange- 
wandt werden, beispielsweise ein Mean-, Median- oder ARMA- 
Filter . 

Des Weiteren ist aus US-Patent 4,463,375 ein Verfahren der 
medizinischen Bildverarbeitung bekannt. Dieses Verfahren gibt 
eine im Rauschen reduzierte Version eines ersten bearbeiteten 
Bildes an, das aus einer Vielzahl von Messungen eines Viel- 
aufnahmensystems, z.B. mittels eines Computertomographiege- 
rats, gewonnen wurde . Im Verfahren wird zuerst aus der Viel- 
zahl von Messungen ein zweites Bild mit hohem SNR erzeugt. 
Dann wird das erste Bild mit einem das Rauschen reduzierenden 
Filter bearbeitet. Das zweite Bild wird mit einem diesem Fil- 
ter komplementaren Filter bearbeitet. Durch die gewichtete 
Kombination der beiden Bildern erhalt man die im Rauschen re- 
duzierte Version des ersten bearbeiteten Bildes. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zum 
Bearbeiten von Magnetresonanzrohdaten sowie ein Magnetreso- 
nanztomographiegerat und ein Magnetresonanzspektroskopiegerat 
anzugeben, die Aufnahmen oder Spektren mit hoher Kanten- 
scharfe erzeugen, ohne dass dabei das Bildrauschen, d.h. das 
SNR, verschlechtert wird. 

Bezogen auf das Verfahren wird diese Aufgabe gemafi der Erfin- 
dung dadurch gelost, dass zur Nachbearbeitung von Magnetreso- 
nanzrohdaten ein Verfahren mit folgenden Verf ahrensmerkmalen 
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durchgefiihrt wird. In einem ersten Schritt werden die Magnet- 
resonanzrohdaten mit einem ersten Filter gefiltert, in einem 
zweiten Schritt werden diese gefilterten Daten f ouriertrans- 
formiert, in einem dritten Schritt bildet man den Betrag der 
fouriertransf ormierten Daten und erhalt so ein erstes Magnet- 
resonanzsignal. In einem vierten Schritt werden die ungef il- 
terten, ursprunglichen Magnetresonanzrohdaten f ouriertrans- 
formiert, und anschliefiend im funften Schritt folgt eine 
zweite Betragsbildung der ungef ilterten fouriertransf ormier- 
ten Daten. Abschliefiend fuhrt man eine gewichtete Kombination 
der beiden Magnetresonanzsignale durch und erhalt so das 
nachverarbeitete Magnetresonanzsignal . 

In anderen Worten besteht das Verfahren aus der gewichteten 
Kombination zweier Magnetresonanzsignale, wobei die Magnetre- 
sonanzsignale aus den gleichen Magnetresonanzrohdaten einmal 
mit und einmal ohne Filterung gewonnen wurden. 
Durch geeignete Wahl des Filters und der Parameter der ge- 
wichteten Kombination extrahiert man eine hinsichtlich der 
medizinisch gewunschten Darstellung optimierte Auf nahme . Ein 
Vorteil des Verfahrens liegt darin, dass das SNR nach der ge- 
wichteten Kombination der beiden Magnetresonanzsignalen er- 
heblich verbessert ist und dass gleichzeitig ein kantenbeton- 
tes Magnetresonanzsignal entsteht. Im Vergleich mit dem Ver- 
fahren nach US-Patent 4,463,375 ist der Aufwand dieses Ver- 
fahrens erheblich reduziert, da es zum einen nur auf den Roh- 
daten einer einzigen Messung basiert und da es zum anderen 
unabhangig von einem Bild mit hohem SNR ist, das aus mehreren 
Messungen erzeugt wird. 

In einer besonders einfachen Erweiterung des Verfahrens wer- 
den die Magnetresonanzrohdaten vor der Erzeugung des zweiten 
Magnetresonanzsignals mit einem zweiten Filter gefiltert. 
Dies ermoglicht eine Einf lussnahme auf das zweite Magnetreso- 
nanzsignal und erlaubt es somit die gewichtete Kombination zu 
verbessern. Die Verwendung der ungef ilterten Magnetresonanz- 
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rohdaten im ursprunglichen Verfahren entspricht der Anwendung 
des Identitatsf liters und ist im Ablauf schneller. 

Eine gewichtete Kombination im Sinne der Erfindung ist dabei 
5 eine mathematische Kombination der Betrage der beiden Magnet- 
resonanzsignale . Diese haben die gleiche Ortskodierung und 
stellen beispielsweise in einer ein- oder mehrdimensionalen 
Pixelstruktur nach der Fouriertransf ormation das Echosignal 
einer Magnetresonanzmessung dar. Die Kombination, meist ist 
10 es eine Addition, erfolgt dabei raumlich konsistent, d.h. es 
werden gleiche raumliche Bereiche kombiniert. Z.B. findet die 
Kombination Pixel fur Pixel statt. Die Gewichtung erfolgt 
liber einen Gewichtsf aktor . Dieser Gewichtsf aktor hangt fur 
einen raumlichen Bereich des einen der beiden Magnetresonanz- 
15 signale vom Inf ormationsgehalt , z.B. der Intensitat, dessel- 
ben raumlichen Bereichs im anderen Magnetresonanzsignal ab. 
Der Gewichtsf aktor kann beispielsweise linear vom Informati- 
onsgehalt abhangen oder jede andere giinstige nichtlineare Ab- 
hangigkeit vom Inf ormationsgehalt des raumlichen Bereichs ha- 
20 ben. Der Gewichts f aktor kann beispielsweise durch eine mathe- 
matische Stuf enfunktion realisiert sein, die unterhalb eines 
vorher festgelegten Schwellwertes des Betragswerts des 
betreffenden Magnetresonanzsignals den Gewichtsf aktor auf ei- 
nen kleinen Wert oder Null setzt. 



30 



Ein Vorteil der gewichteten Kombination liegt darin, dass das 
SNR nach der gewichteten Kombination der beiden Filter erheb- 
lich verbessert ist, und dass durch diese Datenwichtung ein 
kantenbetontes Profil im Magnetresonanzsignal entsteht. 



Das Verfahren bezieht sich auf die Nachverarbeitung von Mag- 
netresonanzrohdaten einer einzelnen Messung und hangt nicht 
davon ab, ob die Rohdaten zu einer Einzelauf nahme gehoren 
oder ob sie zu einer Aufnahme aus einer Auf nahmenserie sind, 
35 die beispielsweise eine zeitliche Entwicklung wiedergeben. 



200209304 



5 

Die Magnetresonanzrohdaten konnen Rohdaten einer Messung ei- 
nes Magnetresonanzspektroskopiegerates sein oder auch Rohda- 
ten einer Aufnahme eines Magnetresonanztomographiegerates 
sein . 

5 

Entsprechend wird die auf ein Magnetresonanztomographiegerat 
bezogene Aufgabe gemaft der Erfindung durch ein Magnetreso- 
nanztomographiegerat gelost, das an ein Verfahren zur Nachbe- 
arbeitung von Magnetresonanzrohdaten mit den obigen Verfah- 
10 rensmerkmalen angepasst ist. Dabei werden die Magnetresonanz- 
rohdaten z.B. als Spinechosignal oder Echosignal bezeichnet. 
^ Die Magnetresonanzsignale entsprechen beispielsweise den ein 
oder mehrdimensionalen Magnetresonanzauf nahmen . 

15 Die auf ein Magnetresonanzspektroskopiegerat bezogene Aufgabe 
wird gemafi der Erfindung durch ein Magnetresonanzspektrosko- 
piegerat gelost, das an ein Verfahren zur Nachbearbeitung von 
Magnetresonanzrohdaten mit den obigen Ver f ahrensmerkmalen an- 
gepasst ist. Dabei werden spektroskopischen Daten, z.B. das 

20 FID-Signal (FID Free Induction Decay), nachbearbeitet . 

In einer besonders gunstigen Ausf iihrungsf orm des Verfahrens 
findet die Anwendung der Filter auf die Magnetresonanzrohda- 
ten nach deren Demodulation statt. Dies erlaubt ein Anwenden 
5 der Filter in ihrer iiblichen Darstellungsf orm im Frequenz- 

raum. Eine Anwendung der Filter nach der Demodulation hat den 
Vorteil, dass die Filter sehr schmalbandig sein konnen. 

In einer besonders vorteilhaf ten Ausf iihrungsf orm ist der 
30 erste Filter ein Tiefpassf ilter , der schnelle Anderungen aus 
den Rohdaten herausf iltert und somit hochf requente Oszillati- 
onen unterdriickt. Beispielsweise kann der Tiefpassf ilter eine 
Art Henning-Filter sein. 




35 



In einer weiteren Ausf iihrungsf orm ist der zweite Filter ein 
Hochpassf ilter, der nur die schnellen Anderungen im Signal 
durchlasst, die beispielsweise auch durch Rauschen verursacht 
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werden. Hochpassgef ilterte Magnetresonanzsignale zeigen eine 
Uberhohung von schnell ansteigenden Kanten. 

In einer besonders vorteilhaf ten Ausfiihrungsf orm, die eine 
Selbstwichtung der Kombination verwendet, wird die gewichtete 
Kombination f olgendermaflen durchgef uhrt : 



C = A + >1 



B 



V^max J 



B 



[G/.l] 



Dabei sind A und B Magnetresonanzsignale, die beispielsweise 
durch Tiefpassf ilterung (A) bzw. Hochpassf ilterung (B) ge- 
wichtet kombiniert werden. A^ ax ist das Maximum des Magnetre- 
sonanzsignals A. Die beiden Parameter X und k bestimmen die 
Beitragsstarke des Magnetresonanzsignals B zum nachbearbeite- 
ten Magnetresonanzsignal C. Der Quotient aus B und A^x be- 
wirkt eine verstarkte Korrektur der Pixel, die einen hohen 
Betrag im Magnetresonanzsignal B aufweisen. Die Werte von X 
und k miissen anhand ihres Einflusses und ihrer optimalen Wir- 
kung auf das Magnetresonanzsignal empirisch ermittelt werden. 
Dabei liegen die Werte von X und k bevorzugt zwischen 1 und 
3. 



Ein Vorteil des Quotienten, der B auf A max normiert, besteht 
darin, dass ein unkontrollierter Signalbeitrag von B in Regi- 
onen, in denen eigentlich nichts zu korrigieren oder zu 
verbessern ist, vermieden wird. 

Die gewichteten Kombination geeignet gefilterter Magnetreso- 
nanzsignale hat den Vorteil, dass in einer einzelnen Magnet- 
resonanzmessung das Rauschen im resultierenden Magnetreso- 
nanzsignal geringer ist als bei einer linearen Addition oder 
Subtraktion der beiden gefilterten Magnetresonanzsignale. 

In besonderen Ausfuhrungsf ormen kann das Verfahren nach der 
Erfindung auf Magnetresonanzrohdaten angewandt werden, die 
einen ein- oder mehrdimensionalen zu untersuchenden Raum er- 
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fassen. Je nach Anwendung werden in der Magnetresonanztechno- 
logie beispielsweise eindimensionale Schnitte, zweidimensio- 
nale Schnittbilder oder dreidimensionale Volumendarstellungen 
wiedergegeben. Die Magnetresonanzrohdaten beinhalten die Di- 
mension der jeweiligen Aufnahme. 

In einer vorteilhaf ten Ausf ilhrungsf orm ist die wahrend des 
Verfahrens durchgef uhrte Fouriertransf ormation in ihrer Di- 
mensionalitat so gewahlt, dass die Dimension der Magnetreso- 
nanzrohdaten erhalten bleibt. Damit kann das Verfahren z.B. 
auf zweidimensionale Magnetresonanzschnittbilder oder auf 
dreidimensionale Volumendarstellungen angewandt werden. Auch 
bei Magnetresonanzspektren ist eine ein- oder mehrdimensio- 
nale Darstellung, beispielsweise in drei Dimensionen, mog- 
lich. 



In einer besonders vorteilhaf ten Ausf uhrungs form des Verfah- 
rens werden zwei ortskodierte Magnetresonanzsignale derart 
kombiniert, dass der Beitrag des einen Magnetresonanzsignals 
zur gewichteten Kombination durch die Multiplikation dieses 
Magnetresonanzsignals mit einem Gewichtsf aktor gebildet wird, 
wobei der Gewichtsf aktor derart von dem anderen Magnetreso- 
nanzsignal abhangt, dass er bei einem groflen Magnetresonanz- 
signal grofier ist als bei einem kleinen Magnetresonanzsignal . 

In einer speziellen Ausf uhrungs form hangt der Gewichtsf aktor 
nichtlinear vom Betrag eines Magnetresonanzsignals ab. Min- 
destens in einer Dimension der Magnetresonanzrohdaten wird 
die gleiche nichtlineare Abhangigkeit verwendet, urn uber alle 
Pixel dieser Dimension den gleichen Effekt beispielsweise auf 
die Bildscharfe zu haben. 

In einer besonderen Ausf uhrungs form des Verfahrens werden die 
Magnetresonanzrohdaten mit mehreren Filtern bearbeitet, bevor 
die verschieden gefilterten Magnetresonanzsignale anschlie- 
fiend mittels gewichteter Kombination zu einem einzigen nach- 
verarbeiteten Magnetresonanzsignal zusammengef asst werden. 
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Mehrere Ausfuhrungsbeispiele des Verfahrens, des Magnetreso- 
nanztomographiegerates und des Magnetresonanzspektroskopiege- 
rates nach der Erfindung werden nachfolgend anhand der Fig. 1 
bis 8 naher erlautert. Es zeigen: 

5 

Fig. 1 ein Flussdiagrainm eines Ausf iihrungsbeispiels des Ver- 
fahrens, das in ein Magnetresonanztomographiegerat 
implementiert wurde, 



* 



10 Fig. 2 ein Flussdiagramm eines Ausf iihrungsbeispiels des Ver- 
fahrens, bei dem zwei gefilterte Magnetresonanzsig- 
nale gewichtet kombiniert werden, 



Fig. 3 eine numerische Simulation der Bearbeitung von Mag- 
15 netresonanzrohdaten anhand eines Spinechosignals, das 

durch ein Rechteckob j ekt im Bildgradienten eines Mag- 
netresonanztomographiegerats hervorgeruf en wird, 



20 



Fig. 4 eine numerische Simulation der Bearbeitung, fur den 
Fall, dass ein Tiefpassf ilter auf das Spinechosignal 
aus Figur 3 angewandt wird, 



Fig. 5 eine numerische Simulation der Bearbeitung, fur den 
Fall das ein Hochpassf ilter auf das Spinechosignal 
\J&5 aus Figur 3 angewandt wird, 

Fig. 6 ein Abszissenausschnitt der Figuren 3 bis 5 zur Ver- 
deutlichung der Wirkungsweise des Verfahrens anhand 
des kantenbetonten Profils, das durch die gewichtete 
30 Kombination erreicht wird, 

Fig. 7 eine VergroJierung der Figur 6 im Bereich einer Kante 
des Rechteckobj ekts, 



35 Fig. 8 



eine numerische Simulation der Bearbeitung von Mag- 
netresonanzrohdaten fur den Fall der Magnetresonanz- 
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spektroskopie anhand eines FID-Signals, das durch die 
Summe dreier Lorentz-Linien gebildet wurde, und 

Fig. 9 Vergleich des durch gewichtete Kombination erhaltenen 
Magnetresonanzspektrums mit dem tiefpassgef ilterten 
Magnetresonanzspektrum, wobei beide Magnetresonanz- 
spektren auf das gleiche Rauschen normiert wurden. 

Die Figur 1 zeigt schematisch den Ablauf des Verfahrens und 
die erforderlichen Komponenten, die innerhalb eines Magnetre- 
sonanztomographiegerates 1 fur den Ablauf des Verfahrens er- 
forderlich sind. In einem Magnetresonanztomographiegerat 1 
mit einem herkommlichen, nicht explizit dargestellten und der 
Ortskodierung dienenden Magnet- und Gradientensystem, werden 
mit einem Hochf requenzsender 3 Hochf requenzimpulse in ein Un- 
tersuchungsobjekt 5 eingestrahlt . Die elektromagnetischen 
Signale, die vom Untersuchungsobj ekt 5 emittiert werden, wer- 
den mit einem Hochf requenzempf anger 7 empfangen. Die Aus- 
gangssignale des Hochf requenzempf angers 7 bilden in Form ei- 
nes Spinechosignals 9 die Magnetresonanzrohdaten . Sie werden 
innerhalb einer Datenverarbeitungsanlage 11 nachbearbeitet . 
Dort werden sie entweder ungefiltert oder mit verschiedenen 
Filtern 13 gefiltert in einer Gewichtungseinheit 15 miteinan- 
der gewichtet kombiniert. Das Ergebnis ist ein nachverarbei- 
tetes Magnetresonanzsignal 17 in Form einer Magnetresonanz- 
auf nahme . 

Die Figur 2 zeigt ein Flussdiagramm eines Ausf uhrungsbei- 
spiels des Verfahrens, bei dem zwei gefilterte Magnetreso- 
nanzsignale in der Datenverarbeitungsanlage 11 gewichtet kom- 
biniert werden. Dabei wird das Spinechosignal 9 erst in einer 
Demodulationseinheit 19 demoduliert, so dass es im k-Raum mit 
der Dimension der Abtastung vorliegt. Anschlieftend wird das 
Spinechosignal 9 zwei Datenverarbeitungseinheiten 21 mit je 
einem Filter zugefiihrt. Dort wird das Spinechosignal das eine 
Mai mittels der Tiefpassf ilterrecheneinheit 23 und das andere 
Mai mittels der Hochpassf ilterrecheneinheit 25 gefiltert. 
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Nach den sich anschliefienden Fouriertransf ormationen 27 und 
den Betragsbildungen 29 erhalt man einerseits die tiefpassge- 
filterte Magnetresonanzauf nahme 31 bzw. die hochpassgef il- 
terte Magnetresonanzauf nahme 33. Diese beiden Magnetresonanz- 
auf nahmen werden in der Gewichtungseinheit 15 durch gewich- 
tete Addition zu einer einzelnen nachverarbeiteten Magnetre- 
sonanzauf nahme 35 kombiniert. Die gewichtete Addition erfolgt 
unter Eingabe von Parametern 37 mittels einer Eingabeeinheit 
39. Die Parameter 37 bestimmen die Gewichtung des Pixels 41 
der hochpassgef ilterten Magnetresonanzauf nahme 33 in Abhan- 
gigkeit vom gleichen Pixel 43 der tiefpassgef ilterten Magnet- 
resonanzauf nahme 31. Bei geeigneter Wahl der gewichteten Ad- 
dition und seiner Parameter weist die nachverarbeitete Mag- 
netresonanzauf nahme 35 sowohl ein hohes SNR als auch eine 
grolie Kantenscharf e auf. 

Im Fall der Verwendung des Verfahrens mit einem Magnetreso- 
nanzspektroskopiegerat ist das Flussdiagramm mit dem der Fi- 
guren 1 und 2 vergleichbar . Nur werden im Unterschied Rohda- 
ten von Magnetresonanzspektren mittels des Verfahrens gefil- 
tert und gewichtet addiert, urn eine hohere Auflosung und ein 
groBeres SNR zu erreichen. 

Die Figuren 3 bis 5 zeigen die Ergebnisse einer numerischen 
Simulation eines Spinechosignals eines Rechteckobj ekts im 
Feldgradienten des Magnetresonanztomographiegerates 1. 
Die Simulation bezieht sich auf einen Schnitt durch ein 
Rechteckobj ekt und verdeutlicht die Auswirkung des Verfahrens 
auf die Kantenscharf e des Objekts. Dieses befindet sich im 
Bereich der 128 Pixel der Simulation. Die Figuren zeigen je- 
weils den Real- und Imaginarteil 51 bzw. 53 des Echosignals 
nachdem ein Filter auf das Echosignal angewendet wurde . Auf- 
getragen ist die Filterf unktion und das komplex-wertige Spin- 
echosignal uber einer k x -Achse im k-Raum. Verschiedene Fil- 
terf unktionen verdeutlichen die Wirkungsweise der Filter auf 
das Echosignal im K-Raum. Des weiteren zeigen die Figuren 3 
bis 7 das Ergebnis der Durchfuhrung einer Fouriertransf orma- 
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tion des gefilterten komplex-wertigen Echosignals mit an- 
schliefiender Betragsbildung . Man erhalt so ein simuliertes 
Messergebnis eines Magnetresonanztomographiegerats, in diesem 
Fall einen eindimensionaler Schnitt durch das Rechteckob j ekt . 
5 Das Schnittbild des Messergebnisses wird uber der entspre- 
chenden raumlichen Achse aufgetragen. 

Die Figur 3 zeigt den Realteil 51 und Imaginarteil 53 des un- 
gefilterten Echosignals sowie das Ergebnis der Fouriertrans- 

10 formation und der Betragsbildung des ungef ilterten Echosig- 
nals, hier als Identitatssimulationsergebnis 57 bezeichnet. 

. In der Simulation wurde dabei kein Filter 55 auf das Echosig- 
< nal angewandt. Man erkennt deutlich die Oszillationen im 

Identitatssimulationsergebnis 57 im Bereich der ansteigenden 

15 und abfallenden Flanken. Diese Oszillationen begrenzen die 
Kantenscharfe in der simulierten (nicht nachbearbeiteten) 
Messung. 

In Figur 4 ist der Einfluss des Tief passrohdatenf ilters 59 
20 auf das Spinechosignal dargestellt. Der Tiefpassrohdatenf il- 
ter 59 wurde in Form eines Henningf ilters F H auf das Spin- 
echosignal angewandt, wobei der Henningf ilter F H 
im eindimensionalen k-Raum folgende Form aufweist: 

^ 5 Ft,h =Fh =1/2x \l + QQs{27u{j-Nl2)lN)] [G/.2] 

Dabei ist N die Anzahl der Abtastpunkte der Aufnahme in einer 
Dimension. Die Filtermodulation folgt einer Cosinus-Funktion, 
deren Maximum urn N/2, also in die Mitte des k-Raums, verscho- 

30 ben ist. Der Index j beschreibt den j-ten Abtastpunkt. Der 

Filter unterdruckt im Echosignal den Randbereich des k-Raums 
mit niedrigerer Intensitat. Dies ist deutlich an der redu- 
zierten Oszillationsamplitude des gefilterten Spinechosignals 
im Randbereich zu erkennen. Das Tiefpasssimulationsergebnis 

35 61 in Figur 4 zeigt eine abgerundete Kante des Rechteckob jek- 
tes. Dies ist auch in den Figuren 5 und 6 zu erkennen. 
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In Figur 5 wird die Simulation flir einen Hochpassrohdatenf il- 
ter 63 durchgef tihrt, wobei der Hochpassf ilter im k-Raum ein 

modulierter, mit einem Offset versehener Henning-Filter F H 
ist : 

5 

^*;h=(f h + Offset)* cos 11 (n{j-N 12)1 N) [G/.3] 

Dabei erfolgt die Modulation des Henning-Filters F H mit ei- 
ner Cosinus-Funktion zur n-ten Potenz, wobei die Minima der 

10 Cosinus-Funktion an den beiden Enden und in der Mitte des k- 
Raums liegen. Die zusatzlichen Minima an den Enden haben den 
^ Vorteil, dass sie keine zusatzlichen Artefakte in der Fou- 

riertransf ormation verursachen. So wird zum Beispiel der so- 
genannte „truncation Artefact* der Fouriertransf ormation un- 

15 terdriickt. 

Dabei hat der Offset-Parameter des modulierten Henningf ilters 
den Wert 0,2 und die Potenz n der Cosinus-Modulation hat den 
Wert 1. Der Hochpassrohdatenf ilter 63 verstarkt die Flugel 
20 des Echosignals im k-Raum. Im Realteil 51 und im Imaginarteil 
53 des gefilterten Echosignals erkennt man die verstarkte 
Amplitude der Echosignaloszillationen in den Randbereichen 
und eine Unterdrtickung im Zentralbereich . Das Hochpasssimula- 
tionsergebnis 65 zeigt ein deutliches Uberschwingen im Be- 
\JK}5 reich der Kanten des Rechteckob j ektes . 

Die Figur 6 zeigt einen Ausschnitt aus den Simulationsergeb- 
nissen im Bereich der Pixel 35 bis 55, d.h. im Bereich der 
ansteigenden Flanke* Man erkennt deutlich die Oszillationen 

30 des Identitatssimulationsergebnisses 57. Die Kante des Recht- 
eckobjekts ist im Falle des Tiefpasssimulationsergebnisses 61 
abgerundet. Das Hochpasssimulationsergebnis 65 zeigt gerade 
im Bereich des Anstiegs eine starke Signaluberhohung . Mit 
Hilfe der gewichteten Addition der beiden gefilterten Simula- 

35 tionsergebnisse kann nun ein kantenbetontes Nachbearbeitungs- 
ergebnis 67 erreicht werden. Dazu wurde eine gewichtete Kom- 
bination nach vorstehend schon angegebener Gleichung 1 durch- 
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gefuhrt. Dabei entsprechen das Tiefpasssimulationsergebnisses 
61 dem Magnetresonanzsignal A und das Hochpasssimulationser- 
gebnis 65 dem Magnetresonanzsignal B. 

Die Parameter X und k werden so eingestellt, dass im Kanten- 
5 bereich eine maximale Scharfe eintritt. Dabei liegen X und k 
meist zwischen dem Wert 1 und 3. Wie stark die Gewichtung des 
Beitrags von B ausfallt, hangt vom Verhaltnis des Betrags im 
jeweiligen Pixel in B zum maximalen Wert A^x der Pixel in A 
ab, d.h. je hoher der Wert des Hochpasssimulationsergebnisses 

10 65 in einem Pixel ist, desto starker ist sein Beitrag. X 

fuhrt dann eine Art lineare Gewichtung durch, wahrend k als 
^ Potenz des Verhaltnisses von B und A den Einfluss dieses Quo- 
tienten gewichtet, also angibt, wie stark die Korrektur von A 
durch das Grofienverhaltnis von B zu A bestimmt wird. Die Pa- 

15 rameter X und k sind in den Figuren 6 und 7 auf die Werte 0, 6 
bzw. 2,3 eingestellt. Neben der Kantenscharfe wird auch das 
SNR des Nachbearbeitungsergebnisses 67 erhoht, das nun mit 
dem SNR des Tiefpasssimulationsergebnisses 61 vergleichbar 
ist. Dies erkennt man an der Unterdriickung des Signals vor 

20 der ansteigenden Flanke, d.h. es liegt ein geringerer Sockel- 
beitrag aufierhalb des Objektes vor. 

Die Figur 7 zeigt einen nochmals vergrofierten Ausschnitt der 
simulierten Messergebnisse . Die Darstellung fokussiert auf 
sjgp den Kantenbereich des Rechteckob j ekts im Bereich der Pixel 42 
bis 55 mit Signalhohen zwischen 10 bis 13. Deutlich erkennt 
man, wie das stark oszillierende Identitatssimulationsergeb- 
nis 57 mittels des Tiefpassf ilters von den hochf requenten Os- 
zillationen befreit wird. Anschlieliend wird durch die gewich- 

30 teten Addition der Eckbereich des Rechteckob j ekts aufgrund 

der hohen Signalbetrage des Hochpasssimulationsergebnisses 65 
aufgefullt. Die Oszillationen des Identitatssimulationsergeb- 
nisses 57 werden auch als „ringing* bezeichnet. Dieser Arte- 
fakt, auch als truncation-Artefact bezeichnet, wird mit Hilfe 

35 des Verfahrens nach der Erfindung unterdruckt. 



2002093*04 



14 

Die Figur 8 zeigt den Ablauf des Verfahrens angewandt auf 
Magnetresonanzrohdaten einer simulierten Magnetresonanzspekt- 
roskopiemessung, in der ein FID-Signal 71 gemessen wird, das 
durch die Summe dreier Lorent z-Linien (73, 75, 77) gebildet 
wird. In Figur 8 sind die Realteile der im Verfahren betei- 
ligten FFT-Spektren dargestellt. In der Simulation wurde ein 
Rauschbeitrag auf das synthetisierte FID-Signal 71 addiert. 
Eine Tief passf ilterung des FID-Signals 71 ergibt das tief- 
passgef ilterte Spektrum 79 und eine Hochpassf ilterung das 
hochpassgef ilterte Spektrum 81. Nach der gewichtete Kombina- 
tion der beiden gefilterten Spektren erhalt man das nachver- 
arbeitete Magnetresonanzspektrum 83. Man erkennt ein deutlich 
reduziertes Rauschen sowie eine hohe Auf losbarkeit der drei 
Lorentz-Linien aufgrund der Anwendung des Verfahrens nach der 
Erf indung. 

In Figur 9 wird das nachverarbeitete Magnetresonanzspektrum 
83 mit dem tief passgef ilterten Spektrum 79 verglichen. Urn den 
Einfluss des Verfahrens zu verdeutlichen, wurden dabei die 
beiden Spektren auf das gleiche Rauschen normiert. Deutlich 
erkennt man die schmaleren Linienbreiten und eine erhebliche 
Steigerung des Kontrasts im nachverarbeitete Magnetresonanz- 
spektrum 83. 
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Patentanspriiche 

1. Verfahren zur Nachverarbeitung von Magnetresonanzrohdaten 
mit folgenden Verf ahrensmerkmalen : 

5 

- Filtern der Magnetresonanzrohdaten mit einem ersten Filter 
(13) , 

- Fouriertransformation (27) der gefilterten Daten, 

10 

- Betragsbildung (29) der f ouriertransf ormierten Daten ergibt 
^ ein erstes Magnetresonanzsignal, 

- Fouriertransformation (27) der Magnetresonanzrohdaten, 

15 

- Betragsbildung (29) der f ouriertransf ormierten der Magnet- 
resonanzrohdaten ergibt ein zweites Magnetresonanzsignal, 

- Bilden eines nachverarbeiteten Magnetresonanzsignals (17) 
20 mittels einer gewichteten Kombination der beiden Magnetre- 

sonanzsignale . 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass die 
^5 Magnetresonanzrohdaten zur Erzeugung des zweiten Magnetreso- 
nanzsignals vor der Fouriertransformation (27) mit einem 
zweiten Filter (13) gefiltert werden. 

3. Verfahren nach Anspruch 1 oder 2, 

30 dadurch gekennzeichnet, dass die 

Magnetresonanzrohdaten nach der Demodulation gefiltert wer- 
den . 

4. Verfahren nach einem der Anspruche 1 bis 3, 

35 dadurch gekennzeichnet, dass die 

Magnetresonanzrohdaten mit einem Magnetresonanzspektroskopie- 
gerat gewonnen werden. 
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5. Verfahren nach einem der Anspriiche 1 bis 4, 
dadurch gekennzeichnet, dass die Mag- 
netresonanzrohdaten mit einem Magnetresonanztomographiegerats 
(1) gewonnen werden. 

6. Verfahren nach einem der Anspriiche 1 bis 5, 

dadurch gekennzeichnet, dass einer der 
Filter ein Tiefpassf ilter ist. 

7. Verfahren nach Anspruch 6, 

dadurch gekennzeichnet, dass der Tief- 
passfilter eine Art Henningf ilter ist. 

8. Verfahren nach einem der Anspriiche 1 bis 7, 

dadurch gekennzeichnet, dass einer der 
Filter ein Hochpassf ilter ist. 

9. Verfahren nach Anspruch 8, 

dadurch gekennzeichnet, dass der Hoch- 
passf ilter eine Art inverser Henningf ilter ist. 

10. Verfahren nach einem der Anspriiche 1 bis 9, 
dadurch gekennzeichnet, dass die 
Magnetresonanzrohdaten Daten eines ein- oder mehrdimensiona- 
len zu untersuchenden Raums sind. 

11. Verfahren nach einem der Anspriiche 1 bis 10, 
dadurch gekennzeichnet, dass die 
Fouriertransformation die Dimensionality der Magnetresonanz- 
rohdaten erhalt. 

12. Verfahren nach einem der Anspriiche 1 bis 11, 
dadurch gekennzeichnet, dass bei der 
gewichteten Kombination die Betrage der Magnetresonanzsignale 
addiert werden. 



13. Verfahren nach einem der Anspriiche 1 bis 12, 
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d a d u r c h 



gekennzeichnet 



dass die 



beiden Magnetresonanzsignale derart gewichtetet kombiniert 
werden, dass der Beitrag eines der Magnetresonanzsignale zur 
gewichteten Kombination durch die Multiplikation dieses Mag- 
5 netresonanzsignals mit einem Gewichtsf aktor gebildet wird, 
wobei der Gewichtsf aktor von dem anderen der Magnetresonanz- 
signale derart abhangt, dass er bei einem groiien Betrag die- 
ses Magnetresonanzsignals grofier ist als bei einem kleinem 
Betrag . 

10 

14. Verfahren nach einem der Ansprtiche 1 bis 13, 
dadurch gekennzeichnet, dass bei der 

™ gewichteten Kombination der Beitrag eines der Magnetresonanz- 
signale nichtlinear vom Betrag des anderen Magnetresonanzsig- 

15 nals abhangt, wobei insbesondere mindestens in einer Dimen- 
sion der Magnetresonanzrohdaten die gleiche nichtlineare Ab- 
hangigkeit verwendet wird. 

15. Verfahren nach einem der Ansprtiche 1 bis 14, 

20 dadurch gekennzeichnet, dass die ge- 
wichtete Kombination der beiden Magnetresonanzsignale durch 
zwei Parameter A, und k in folgender Weise das nachverarbei- 
tete Magnetresonanzsignal ergibt: 




C = A + A B , wobei 



A eines der Magnetresonanzsignale ist, 

B das andere der Magnetresonanzsignale ist, 

Amax das Maximum des Magnetresonanzsignals A ist und 

C das nachverarbeitete Magnetresonanzsignal ist. 



30 



16. Verfahren nach einem der Ansprtiche 1 bis 15, 
dadurch gekennzeichnet, 
rameter X im Bereich zwischen 0 und 3 liegt. 



dass der Pa- 



35 



17. Verfahren nach Anspruch 15 oder 16, 
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dadurch gekennzeichnet, dass der Pa- 
rameter k im Bereich zwischen 0 und 3 liegt. 

18. Verfahren nach einem der Anspriiche 1 bis 17, 
dadurch gekennzeichnet, dass mehr 
als zwei Magnetresonanzsignale aus den Magnetresonanzrohdaten 
mittels Filterungen gewonnen und gewichtet zu einem nachver- 
arbeiteten Magnetresonanzsignal kombiniert werden, 

19. Magnetresonanztomographie-Gerat (1), das an ein Verfahren 
zur Nachverarbeitung von Magnetresonanzrohdaten nach einem 
der Anspriiche 1 bis 18 angepasst ist. 

20. Magnetresonanzspektroskopie-Gerat , das an ein Verfahren 
zur Nachverarbeitung von Magnetresonanzrohdaten nach einem 
der Anspriiche 1 bis 18 angepasst ist. 
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Zusammenf assung 

Verfahren zum Bearbeiten von Magnetresonanzrohdaten, Magnet- 
resonanztomographiegerat und Magnetresonanzspektroskopiegerat 

Bei einem Verfahren zum Erzeugen einer Abbildung aus Magnet- 
resonanzrohdaten eines Magnetresonanztomographiegerats (1) 
oder eines Spektrums aus Magnetresonanzrohdaten eines Magnet- 
resonanzspektroskopiegerats werden die Rohdaten mit einem 
Hochpass- und mit einem Tiefpassf ilter gefiltert. Anschlie- 
fiend wird aus den beiden Betragen der gefilterten, fou- 
riertransformierten Magnetresonanzrohdaten die Abbildung 
oder das Spektrum mittels gewichteter Kombination erzeugt. 
Das Verfahren zeichnet sich dadurch aus, dass das Signal-zu- 
Rausch-Verhaltnis zusammen mit einer Erhohung der Kanten- 
scharfe verbessert werden kann. 



Fig 1 
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FIG 3 
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FIG 8 
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FIG 9 




